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The reliable extraction of information from the two- proton 
correlation functions measured in heavy-ion reactions is a 
long-standing problem. Recently introduced imaging tech- 
niques give one the ability to reconstruct source functions 
from the correlation data in a model independent way. We 
explore the applicability of two-proton imaging to realistic 
sources with varying degrees of transverse space-momentum 
correlations. By fixing the freeze-out spatial distribution, we 
find that both the proton images and the two-particle corre- 
lation functions are very sensitive to these correlations. We 
show that one can reliably reconstruct the source functions 
from the two-proton correlation functions, regardless of the 
degree of the space-momentum correlations. 



The sensitivity of the two-proton correlations to the 
space-time extent of nuclear reactions was first pointed 
out by Koonin |l| and later emphasized by many au- 
thors [0-k| • Since then, measurements of the two-proton 
correlations have been used along with pion HBT data 
as a probe of the space-time properties of the heavy- 
ion collisions (for the review of recent experimental re- 
sults of two-particle interferometry see [g|-|8| and refer- 
ences therein). A prominent "dip+peak" structure in 
the proton correlation function is due to the interplay of 
the strong and Coulomb interactions along with effects of 
quantum statistics. Because of the complex nature of the 
two-proton final state interactions only model-dependent 
and/or qualitative statements were possible in proton 
correlation analysis. Typically [B-13|, the proton source 
is assumed to be a chaotic source with gaussian profile 
that emits protons instantaneously. For simple static 
chaotic sources, it has been shown pl-pj that the height of 
the correlation peak approximately scales inversely with 
the source volume. Heavy-ion collisions are complicated 
dynamic systems with strong space-momentum correla- 
tions (such as flow) and a nonzero lifetime. Hence the va- 
lidity of the assumptions behind such simplistic sources 
is questionable. In order to address the limitations of 
this type of analysis (and to incorporate collective ef- 
fects), some authors ||[[3],[L4| utilize transport models to 
interpret the proton correlation functions. Although this 



approach is a step in a right direction, it is still highly 
model-dependent. 

Recently, it was shown that one can perform model- 
independent extractions of the entire source function 
S(r) (the probability density for emitting protons a dis- 
tance r apart) from two-particle correlations, not just 
its radii, using imaging techniques flSTLq] . Furthermore, 
one can do this even with the relatively complicated pro- 
ton final-state interactions and without making any a pri- 
ori assumptions about the source geometry or lifetime, 
etc. First results from the application of imaging to the 
proton correlation data can be found in refs. jl^ , |l(||l9| . 
While these results look promising, tests of the imaging 
technique have only been performed on static (Gaussian 
and non-Gaussian) sources. It is important to under- 
stand the limitations and robustness of this technique 
especially in the light of the ongoing experimental pro- 
gram at SIS, AGS and SPS as well as upcoming experi- 
ments at RHIC. 

In this letter, we will study the applicability of pro- 
ton imaging to realistic sources with transverse space- 
momentum correlations. In particular we explore how 
?T — pt correlations directly affect the proton sources 
and, hence, the shapes of the experimentally observable 
correlation functions. Here ?t and pt are the transverse 
radius and transverse momentum vectors respectively of 
a proton at the time when it decouples from the system 
(freeze-out). It has been argued in the pion [E0[ and pro- 
ton p) cases that the apparent source size (or the effective 
volume) decreases as collective motion increases. We will 
verify this expectation and show that one can reliably re- 
construct the source function, even in the presence of ex- 
treme space-momentum correlations. The outline of this 
letter is as follows. First, we briefly describe the imaging 
procedure used to extract the source function from ex- 
perimental correlations. We will discuss proton sources 
but most of our arguments and conclusions are valid for 
any two-particle correlations. Next, we describe how we 
implement the varying degrees of space-momentum cor- 
relations using the RQMD model. Finally, we will discuss 
the influence of these correlations on the proton correla- 
tion functions and imaged sources. Since we can also 
construct the sources directly within RQMD, this serves 
as a more demanding test of the imaging procedure than 



has been performed to date. 

With imaging, one extracts the entire source func- 
tion S(r) from the two-proton correlation function, C(q). 
Here the source function is the probability density for 
emitting protons with a certain relative separation in the 
pair Center of Mass (CM) frame. The source function 
and the correlation function are related by the equation 



C{q) - 1 = 4tt / drr 2 K(q, r)S(r) 



(1) 



In eq. (pi), q = \\f{vi ~ P2) 2 is the invariant relative 
momentum of the pair, r is the pair CM separation after 
the point of last collision, and K is the kernel. The kernel 
is related to the two-proton relative wavefunction via 



#M = ~^(2j + l) 



) tf.(r) -1 



(2) 
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Here cn a are the relative proton radial wavefunctions 
for orbital angular momenta £,£', total angular momen- 
tum j, and total spin s. In what follows, we calcu- 
late the proton relative wavefunctions by solving the 
Schrodinger equation with the REID93 [g^l nucleon- 
nucleon and Coulomb potentials. 

Because (|lj) is an integral equation with a non-singular 
kernel, it can be inverted |p.5|l6 |. To perform the inver- 
sion, we first discretize eq. (|lj), giving a set of linear equa- 
tions, Ci — 1 = Y) j—x KijSj, with N data points and M 
source points. Given that the data has experimental er- 
ror ACj, one cannot simply invert this matrix equation. 
Instead, we search for the source vector that gives the 
minimum % 2 : 
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£ 



(Ci-l 



2_/j=i KijSj 



(AC,) 2 



(3) 



The source that minimizes this \ 2 is (i n niatrix notation) : 

S ={K T BK)- 1 K T B{C -I) (4) 

where K T is the transpose of the kernel matrix and 
B is the inverse covariance matrix of the data, B^ = 
Sij j (ACi) 2 . In general, inverse problems such as this one 
are ill-posed problems. In practical terms, small fluctu- 
ations in the data can lead to large fluctuations in the 
imaged source. One can avoid this problem by using the 
method of Optimized Discretization discussed in refer- 



ence 1 16 . In short, the Optimized Discretization method 
varies the size of the r-bins of the source (or equivalently 
the resolution of the kernel) to minimize the relative error 
of the source. 

The source function that one reconstructs is directly 
related to the space-time development of the heavy-ion 
reaction in the Koonin-Pratt formalism |qLbiIP3| : 



S(r,q)= j dVl r I dhdh [ d 3 R 

J Ait J J 



Un J J (5) 

xD{R + f/2, h-q)D{R - f/2, t 2 : -q), 

where q = h(pi — P2), making q — \q\. Here the D's 
are the normalized single particle sources in the pair CM 
frame and they have the conventional interpretation as 
the normalized phase-space distribution of protons after 
the last collision (freeze-out) in a transport model. In 
computing S(r, q) in a transport model, one does not 
need to consider the contribution of large relative mo- 
mentum {q > <7cut) pairs to the source as the kernel cuts 
off the contribution from these pairs. The kernel does 
this because it is highly oscillatory while the source varies 
weakly on the scale of these oscillations and the integral 
in (pi) averages to zero. We can estimate q cu t directly 
from the correlation function as g cu t is roughly the mo- 
mentum where the correlation goes to one. Nevertheless, 
for the imaging in (fy) to be unique one must require that 
the q dependence of the correlation comes from the kernel 
value alone and eq. (ph seems to indicate that the source 
itself has a q dependence. In practice S(r, q) only has a 
weak q dependence for q ^ g cu t and this dependence may 
be neglected |l|,|. 

Since S(r) is the probability density for finding a pair 
with a separation of emission points r, one can compute 
it directly from the freeze-out phase-space distribution 
given by some model. First one scans through this freeze- 
out density of protons, then histograms the number of 
pairs in relative distance in the CM, and finally normal- 
izes the distribution: An J drr 2 S (r) = 1. As mentioned 
above, only low relative momentum pairs may enter into 
this histogram as the kernel cuts off the contribution from 
pairs with q > q cut . 

In our studies we used the Relativistic Quantum 
Molecular Dynamics (RQMD) model p4| . It is a semi- 
classical microscopic model which includes stochastic 
scattering, classical propagation of the particles. It in- 
cludes baryon and meson resonances, color strings and 
some quantum effects such as Pauli blocking and finite 
particle formation time. This model has been successfully 
used to describe many features of relativistic heavy-ion 
collisions at AGS and SPS energies. Our approach is as 
follows: first we take the freeze-out phase space distribu- 
tions generated by RQMD and we alter the orientation of 
transverse momentum relative to the transverse radius, 
obtaining a subset of the phase space points. Following 
this, we use the Lednicky-Lyuboshitz f|,g5| method to 
construct the proton-proton correlation function. This 
method gives a description of the final state interactions 
between two protons, including antisymmetrization of 
their relative wave function. Finally, using the imaging 
technique described above we compute the proton source 
functions. We used w 4000 simulated events of 4 GeV/A 
Au-Au reactions with impact parameter b < 3 fm. We 
utilized only pairs in the central rapidity region with 



\y\ < 0.3 and applied no cut on transverse momentum 

Pt- 

We consider three different degrees of alignment be- 
tween the transverse position fr and the transverse mo- 
mentum pr of each proton used to construct the correla- 
tion function. These alignments arc implemented in the 



same manner as in reference 26 



1. We orient pt at a random angle with respect to fr- 
We refer to this as the random case. One can think 
of this case as being "thermal" as the transverse 
flow component is completely removed. 

2. We do not change the orientation of pr- We refer 
to this case as the unmodified case. 

3. We align pr with Tt and refer to this as the aligned 
case. One can think of this case as one of extreme 
transverse flow. 

Note that the rotation occurs in the rest frame of the col- 
liding nuclei. In all cases, we only rotate pr so these pro- 
cedures do not change the spatial distribution at freeze- 
out. However, it is clear that these procedures do change 
the phase-space density. 




FIG. 1. Upper panels: Original (solid symbols) and re- 
stored (histogram) proton correlation functions for different 
degrees of space-momentum alignment. The error bars in the 
original correlation functions are from the limited statistics 
of the RQMD runs. Lower panels: Reconstructed sources 
(histogram) and source computed directly in RQMD (solid 
symbols) for the different degrees of alignment. 

Upper panels on Fig. [j\ show the correlation functions 
for the three different cases. It is clear from the Fig- 
ure that the degree of space-momentum correlation has 



a strong influence on the correlation function: the peak 
height of the correlation function changes from about 
1.45 for unmodified RQMD to about 2.2 for the aligned 
case and 1.2 for randomized case. We would like to stress 
again, in all cases the spatial part of the source e.g. "ra- 
dius" or "volume," remains unaltered as does the trans- 
verse momentum spectrum and rapidity distribution of 
protons. Hence, the upper panels of Fig. 1 illustrate the 
danger of ignoring space-momentum correlations when 
analyzing correlation data. 

On lower panels in Fig. [j] we show the proton sources 
obtained with the help of the imaging procedure outlined 
above. Notice that, as the degree of alignment increases 
(going from right to left), that the source function be- 
comes narrower and higher. 



Aligned 



Randomized 





FIG. 2. Sketch of single particle sources at central rapidity, 
looking in the beam direction for the aligned case (left) and 
random case (right). The small circles represent the protons 
and the arrows represent their transverse momentum vectors. 
The areas outlined with the dashed lines represent regions 
where we find pairs with small relative momentum. 

One can understand this shift to lower separations in 
the way sketched in Fig. ||. In the aligned case, it is more 
probable that nearby protons have a small relative mo- 
mentum q. In the random case, any pair can have a small 
q, regardless of their separation. Given that the kernel 
cuts off contributions from pairs with larger q, we expect 
that the aligned case will have a narrower source than 
the unmodified case and the unmodified case will have a 
narrower source than the random case. Also, in Fig. fil 
we have shown the sources constructed directly from the 
RQMD freeze-out distribution following eq. (ph. In these 
sources, we considered all pairs with a relative momen- 
tum smaller than q cut = 60 MeV/c. We explored a range 
of <?cut of 60 to 100 MeV/c, all beyond the point in the 
correlation where it is consistent with one, and found no 
cutoff dependence. 

In all cases we see a general agreement of the im- 
aged sources with the low relative momentum sources 
constructed directly from RQMD. In order to check the 
quality of the imaging and numerical stability of the in- 



version procedure, the two-proton correlation functions 
were calculated using the extracted relative source func- 
tions shown on lower panels in Fig. fy as an input for 
eq. (m) . The result of such "double inversion" procedure 
is shown on upper panels in Fig. nl with solid circles. The 
agreement between the measured and reconstructed cor- 
relation function is quite good, confirming that imaging 
produces numerically stable and unbiased results. 

In conclusion, we have explored the applicability of 
proton imaging to realistic sources with transverse space- 
momentum correlations. By fixing freeze-out spatial dis- 
tribution and varying the degree of transverse space- 
momentum correlation we found that both the images 
and the two-particle correlation functions are very sensi- 
tive to these correlations. In particular, we have shown 
that the source function narrows (i.e. the probability 
of emitting pairs with small relative separation grows) 
and the peak of the proton correlation function in- 
creases as the degree of alignment increases. Finally, 
we have demonstrated that one can reliably reconstruct 
the source functions even with extreme transverse space- 
momentum correlations. We would like to point out that 
the effects of space-momentum correlations should be 
even more pronounced in the shapes of three-dimensional 
proton sources. Note that three dimensional proton 
imaging is now possible J27) . 

An important direction for the future is a detailed 
study of the change of the phase-space density and en- 
tropy (extracted from imaged sources p5[) with the vary- 
ing degree of space-momentum correlation. Such work 
should provide information complementary to ongoing 
studies in the pion sector p8| ]. 
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